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We calculate within a semiclassical approximation the autocorrelation function of cross sections. 
The starting point is the semiclassical expression for the diagonal matrix elements of an operator. 
For general operators with a smooth classical limit the autocorrelation function of such matrix 
elements has two contributions with relative weights determined by classical dynamics. We show 
how the random matrix result can be obtained if the operator approaches a projector onto a single 
initial state. The expressions are verified in calculations for the kicked rotor. 



I. INTRODUCTION 



Quantum systems whose classical limit is chaotic show 
fluctuations in cross sections and eigenvalue positions 
whose statistical properties seem to fall into a few uni- 
versality classes Among the many measures that 
have been applied to characterize these statistical fea- 
tures, much attention has been given to two-point corre- 
lation functions since they can under certain assumptions 
be related to the classical dynamics [Q. For the case of 
spectra of bounded systems this has worked remarkably 
well and in addition one of the main predictions of the 
semiclassical analysis, the existence of long range corre- 
lations due to periodic orbits has been confirmed 
many times [^]J^ . 

More recently investigations of the statistical be- 
haviour of directly observable quantities, such as cross 
sections, have been worked out within the nonlinear 
sigma model for disordered systems The corre- 

lation function was found to have two contributions, a 
Lorentzian and a derivative of a Lorentzian with respect 
to its parameter. The ratio between the two terms is 
fixed and depends on symmetry only. Since there is no 
semiclassical expression for the individual wave functions 
from which the cross sections could be calculated, the 
derivation of such correlation functions within semiclas- 
sics poses a serious challenge. A first step in this direc- 
tion was undertaken by Agam , who exploits quantum 
properties of the matrix elements and does not use previ- 
ously established formulas for diagonal matrix elements 
pO| , pl| . The derivation presented here is similar in spirit 
but starts from the semiclassical expression for diagonal 
matrix elements and specializes to the case of the cross 
section in the end. In particular, we show how the rela- 
tive weight between the two contributions to the correla- 
tion function can be changed. The final expressions are 
compared with data for cross sections in an open kicked 
rotor model. 

In section 2 we present the semiclassical derivation of 
the correlation function between cross sections. This cal- 



culation is actually straightforward and closely patterned 
after calculations for other two point correlations In 
section 3 we discuss the limit that has to be taken in the 
observable to arrive at the correlation function for cross 
sections. In section 4 we discuss numerical simulations 
for an open kicked rotator. Some concluding comments 
are given in section 5. 

II. SEMICLASSICAL CORRELATION 
FUNCTIONS FOR SMOOTH OPERATORS 

Quantum cross sections for the transition from an ini- 
tial state \i) to a set of final states \n) which are eigen- 
states of a Hamiltonian, H\n) = En\n), are proportional 
to 

a{E) cx ^{n\D\i) {i\D\n)6{E - H) (1) 

n 

where D denotes the dipole operator. Using the projec- 
tion operator 

A = D\i){i\D (2) 

and Greens function G = \/{E — H) the cross section 
becomes (dropping the proportionality factors implied in 

(1)) 

aA{E)^--lmiv{GA). (3) 

If the operator A is sufficiently smooth and has a non- 
singular classical limit, there is a semiclassical expression 
for a{E) which naturally divides into two pieces | p^ , pl| , 

af\E)^a^{E)+aAji{E). (4) 

The first term, a smoothly varying background contri- 
bution from paths of 'zero length', is determined by in- 
tegration over the energy shell of the Wigner transform 
Aw of the observable with the measure d/i = d'^pd'^q/h'^ 
(in d-degrees of freedom). 
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ao{E) = J dnAwS{E-H) 



(5) 



The second term describes the fluctuations around it and 
is determined by classical periodic orbits, 



/h 



(6) 



where Sp, Tp and Wp are the action, period and weight 
of the periodic orbit p, respectively, and 



Ap — 



dtAwipit),q{t)) 



(7) 



is the integral of the observable (again the Wigner trans- 
form of A) over the periodic orbit p. This result is de- 
rived under the condition that the observable is suffi- 
ciently smooth so that it does not affect the stationary 
phase evaluation that singles out periodic orbits. Wigner 
transformations of projection operators are critical in this 
respect since they typically approach delta functions in 
momenta as h goes to zero, i.e. become rather singular. 
This is the main limitation that prevents a direct appli- 
cation of the above expression to the calculation of the 
autocorrelation function of the cross sections. We there- 
fore adopt the following strategy: We first calculate the 
auto-correlation function for aAiE) within a semiclassi- 
cal approximation for a smooth observable A and then 
discuss the limit of a singular operator. 

The object we want to calculate is the normalized auto- 
correlation function of the fluctuations around the mean 
cross section (7o{E), 

C{e) = {aAji{E + e/2)aAji{E - el2))/{cT^f (8) 

where (. . .) denotes an average over energy. The energy 
scale is set quantum mechanically by the mean spacing 
A between neighboring levels, calculated from the mean 
density of states 



Pq = 1/A = j dti5{E - H) . 



(9) 



The associated time scale is the Heisenberg time Th — 
2TThpQ — h/A. The scale for the observable is set by the 
phase space average of its Wigner transform. 



A 



J dpAwSjE - H) 
J dp,6{E - H) 



(10) 



so that ctq = Apa = A/A. 

Substituting the fluctuating part from the periodic or- 
bits gives 

A2 

C[e) = ^ 



2A^^h^ 

Re/^A;Ap,u;>p.e^(^''(^+^/2)-V(^^--/2))A^ 
\p,p' / 



The classical action in the exponent can be expanded for 
small e, 

Sp{E + e/2)-Sp,{E-el2)^ 

Sp{E) - Sp,{E) + e{Tp + Tp,)/2 . 

In the diagonal approximation ^] the correlation func- 
tion becomes 



Cdi: 



iTpe/h 



(11) 



the factor g accounts, as usual, for degeneracies in case of 
real symmetric Hamiltonians, g — 2, while for hermitian 
Hamiltonians 5 = 1. 

For the next step we use the periodic orbit sum rule 
|l2| and follow the steps in jisj . The periodic orbits pro- 
liferate exponentially with time so that the sum on p can 
be replaced by an integral over time. The density of or- 
bits is given by e'^-^/T with p the topological entropy and 



their weight by |?«p| 



-AT 



with A the Lyapunov ex- 



ponent. The difference between the Lyapunov exponent 
and the topological entropy is the classical escape rate. 



T = X-p. 



(12) 



The integrals Ap of the observable along the orbit vary 
considerably among orbits of similar length. The high 
density of periodic orbits allows to capture this proba- 
bilistically through the distribution P{A) of values ob- 
tained for all orbits with periods in a small interval 
around T. If the correlations in the classical dynamics 
fall off sufficiently rapidly, the distribution will be Gaus- 



PiA) 



l-KSA 



-(A-Af/sl 



(13) 



with a mean A — AT following from ergodicity and a 
variance s\ = oT that increases linearly with time. With 
this distribution function and the assumption that there 
are no correlations between weights Wp and observables 
Ap, the mean square average of the Ap's from orbits with 
periods near T changes with time like 



{\ApY){T)^ A'T' + aT 



(14) 



Thus summing the contributions of orbits in the diagonal 
approximation we obtain 

.gA2 

i-Th (n-X+ie/n)T , 

Re j dT- [AT^ + aT 

For bounded systems {p = A) and for e = we obtain the 
expression for the variance of matrix elements derived in 
Ip^ and tested and verified in many situations [p^-|l6t . If 
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the system is sufficiently open so that TTh ^ 1 the con- 
tributions from the ofF-diagonal terms can be neglected 
and the integration continued to infinity (as explained in 
the next paragraph). Then 



2A n^n^ 



-Re a 



r - is/h 



A 



1 



(r-ie/n)2 

(15) 



With the energy and time scales mentioned above, i.e. 
e = eA and F = 2ttT/Th = TA/h, the correlation func- 
tion becomes 



Cie) 



1 r2 - 



^ \a th + £2 2tt (r2 -1- £2 



(16) 



As in the calculation of Fyodorov and Alhassid the 
correlation function has two terms, a Lorentzian and 
a derivative of a Lorentzian with respect to the width. 
However, in contrast to their formula, where the relative 
weight between the two terms was fixed, it here depends 
on the observable and the Heisenberg time. This point 
will be taken up again in the next section. 

The above derivation is based on the usual assump- 
tions on the diagonal approximation, the validity of a 
periodic orbit sum rule and the replacement of a sum 
over orbits by an integral in time. As a consequence de- 
viations can be expected for short times where isolated 
periodic orbits dominate, an effect that should be partic- 
ularly noticeable near bifurcations. Deviations from the 
diagonal approximations are strong for the orthogonal 
ensemble and absent for the unitary ensemble, at least 
up to the Heisenberg time |jl^. The extension of the 
time integration up to infinity rather than the Heisen- 
berg time is justified if the system is very open, i.e. if 
TTh ^ 1, so that the corrections due to off-diagonal 
terms for times beyond the Heisenberg time can be ne- 
glected. When approaching bounded systems T vanishes 
and the corrections have to be taken into account. The 
final expression is thus reasonable only for sufficiently 
open systems where T/h ^ A. In this limit the contri- 
butions from long orbits are quickly suppressed and the 
differences between the unitary and orthogonal ensem- 
bles should disappear, except of course for the factor of 
2. It is possible to go beyond this by assuming that the 
form factor is the random matrix form factor times an 
exponential damping, as suggested also by Alhassid and 
Fyodorov |,|l|l. 

The derivative of the Lorentzian in the second term, 
weighted by the mean of the operator can be traced back 
to the autocorrelation function of the density of states 
(without operator) as calculated earlier [|l^. The width 
r that enters here is the classical escape rate, since that 
is what determines the modification of the classical sum 
rule. The prediction then is that the quantum resonances 
have half that width, since it is the probabilities and 
not the amplitudes that have to follow the classical be- 
haviour. In many cases, especially with a finite num- 



ber of channels, the situation presumably is more com- 
plex, since the quantum resonances have a distribution of 
widths and it is not clear which quantity (average width, 
maximum of distribution, longest life time etc) dominates 
the form factor |Q . In the semiclassical limit of an infi- 
nite number of channels some information can be drawn 
from the results of Fyodorov and Somniers for random 
matrix models with fixed transmission [ pT| and from the 
distribution of resonance widths calculated by Haake et 
al. 22 1 : in both cases a single width parameter, given by 



half the classical escape rate in the first case and by half 
the gap in the distribution in the second case suffices to 
describe the correlation function. 

Finally, we note that if the classical correlations do 
not decay sufficiently rapidly the average of the square 
of the integrals along the orbits can be expressed as an 
integral over the correlation function as in the previous 
calculation of matrix elements [llSl . 



III. THE LIMIT OF PROJECTION OPERATORS 

The final result for the correlation function given in Eq. 
( p^ ) has the two functional dependencies also identified 
in Alhassid and Fyodorov 0, but the relative weighting 
depends on the observable. Even worse, the first term 
contains Th in the denominator and therefore seems to 
vanish in the semiclassical limit where Th diverges. Then 
the correlation function is of the form of a derivative of 
the Lorentzian only. So how can one obtain Alhassid and 

2 

Fyodorov's result a/ A Th = 1 within this semiclassical 
approach? 

The key to the problem is the observation that the 
semiclassical approximation assumes the observable to 
have a nonsingular classical limit whereas the quantum 
cross sections are obtained from observables that are pro- 
jectors on the initial state (weighted with the dipole op- 
erator). The Wigner transform of a projector is itself 
a function of Planck's constant and becomes singular in 
the semiclassical limit. More specifically, one can visual- 
ize the Wigner transform of a projector as a characteris- 
tic function that in d dimensions lives on a phase space 
cell of volume h'^, since that is the phase space volume 
occupied by a single state. The observables for which 
the semiclassical trace formula was derived were smooth 
with a non-singular limit, that is to say they covered an 
increasing number of quantum states in the semiclassi- 
cal limit. This smearing over many states suppresses the 
first term in the correlation function. 

It is possible to estimate the consequences of this obser- 
vation on the operator A in a simple model of a uniformly 
damped quantum map. Consider a 2-d chaotic map on 
a finite phase space and its quantum representation by 
an iV X unitary operator U . The dimension A'' of U , 
Planck's constant h and the volume JIq of the classical 
phase space are connected by h = flo/N. Damping is 
introduced uniformly everywhere in phase space and on 
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all quantum states. The projector is modelled by an ob- 
servable that takes on the value qq in some part of phase 
space of area Qa and vanishes everywhere else. 

On the classical side the integrals of the observable 
along the periodic orbit are replaced by sums over the 
points of the orbit. If there are no correlations between 
different time steps, the average value of the Ap over 
all orbits with period n is given by (Ap) = naop, where 
p — ^Ia/^o is the probability for a randomly chosen point 
to lie in the phase space area where the observable does 
not vanish. The second moment of the distribution is 
given by 

{Al)^ayn^ + alpil-p)n. (17) 

By comparison with Eq. ( [T^ ) we read off A = aop and 
a = aop(l — p). The Heisenberg time is Th — N, the 
dimension of the Hilbert space. Therefore, 

a _ alp{l-p) 



a'Th 



aw 



I-P 
pN 



(18) 



With an escape rate f expressed in units of the Heisen- 
berg time the correlation function then becomes 



1 



piV r2 -h£2 2tT (r2 -f £2)2 



(19) 



This expression clearly shows the suppression of the first 
term in the semiclassical limit of large TV if p is fixed. 

However, if the observable is a projector onto a single 
state its Wigner transform should localize on a cell of 
phase space volume h. Thus, ^Ia — ^o/N and p — 1/N, 
so that the product pN = 1. Except for the tiny cor- 
rection 1/N to the correct ratio of 1 this is the result 
obtained by Alhassid and Fyodorov 0. In addition, it 
suggests a way to modify the relative weighting between 
the two terms: consider transitions not from a single state 
but from M states in an incoherent superposition, i.e.. 



A = Dj2Wi){m\D . 



(20) 



Then the phase space area covered by A will increase 
M-fold and the weight of the Lorentzian will decrease 
correspondingly. 



1 T'-e' 



n \MT^+£2 27r(r2 + £2)2 



(21) 



(where it is assumed that M <C N). The resolution to 
the problem of the relative weight between the two terms 
in the correlation function posed at the beginning of this 
section thus is that by focusing on initial states which 
are projectors the classical observable also depends on h, 
and this ?i-dependence influences the classical quantities 
as well. Indeed, if the classical observable becomes very 
localized in phase space, it is rarely visited, the return 
time becomes large and the variance increases much more 
slowly, on a timescale also set by the return time. In this 
way the Heisenberg time enters the classical quantities. 



IV. NUMERICAL TESTS ON THE KICKED 
ROTOR 



The model of the previous section clearly stretches the 
applicability of the semiclassical expressions for the ma- 
trix elements to their limits and requires numerical tests. 
We use the kicked rotor in the momentum space quanti- 
zation of Izrailev |23| for this purpose. 



N-l 



Unm - jyYl 



^-iKNV(2-Kl/N) ^-2i-Kl{n-m)/N ^-i2-Km^ /N 



1=0 



(22) 



where [|l8[ V{(j)) = cos — sin(20) is the kicking potential. 
This model is known to belong to the unitary universal- 
ity class {g — 1) due to the second term in the poten- 
tial that breaks the conjugation symmetry ||2^ , ^ . If the 
kicking strength is sufRciently large the correlations de- 
cay very quickly and one ends up with essentially random 
values for the momenta. The calculations were done for 
A' = 7 1^ and matrix sizes N = 101, 201, 401, 801 
and 1601. The initial states were taken to be momentum 
eigenstates of the unperturbed map. From the eigen- 
states {v) and eigenphases </)^ a cross section was formed 
with the operator A from Eq. (EG) according to 



a(0) = Re5]- 



exp(i( 



r/2) 



(23) 



The damping r/2 is uniform for all eigenstates and mod- 
els the coupling of the system to a continuum (for a dis- 
cussion of this point see the previous section and the 
references cited there). The role of the energy is now 
taken over by the phase (f) and the mean and correlation 
function are calculated over the periodicity interval 27r. 
The mean separation between eigenphases is A = 2tt/N. 
Figure 1 shows the correlation function obtained for a dif- 
ferent number of initial states and constant damping. For 
a single initial state the contribution from the derivative 
of a Lorentzian is barely noticeable but as the number of 
initial states increases the deviations from the Lorentzian 
become larger. In particular, the correlation function de- 
velops a zero for M/27r > T. Fitting Eq. (16) with the 
coefficient of the Lorentzian and F as free parameters 
yields a broadening F « F independent of M and the 
coefficient of the Lorentzian in Eq. ( [l6| ) is proportional 
to M^^ as expected from Eq. (]2l|). 

A random matrix calculation |p5|| shows that the quan- 
tity that controls the relative weight between the two 
terms is the ratio between the variance s\ of the ma- 
trix elements and the square of the average for differ- 
ent matrix sizes. This quantity, {sa/{A))'^ should equal 
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(1 — p)/pN. As shown in Figure 2, the renormahzed 
quantity {sa/ {A))'^pN follows the expected 1 — p be- 
haviour rather closely, albeit with large fluctuations. It 
thus seems that the assumptions that enter in our semi- 
classical derivation of the correlation function can be sat- 
isfied in chaotic systems. 



V. FINAL REMARKS 

We have shown how within a semiclassical approxi- 
mation correlation functions for cross sections in open 
systems can be calculated. The calculations could be 
supported with simulations in the standard map and in 
particular the changes in the relative weight between the 
two contributions to the correlation function could be 
demonstrated. 

Actually, the calculation works much better than can 
reasonably be expected: it is well known that the cal- 
culation of wave function usually requires in addition to 
periodic orbits also recurrent orbits ||2^,|2^. Their im- 
portance depends on the width of the initial state [ pT| . 
It may happen, however, that in a statistical sense the 
differences between the contributions from recurrent and 
periodic orbits cancel. A related problem concerns the 
higher moments of the distribution and thus the form 
of the full distribution. Even in the singular limit of 
small p considered in the model, the distribution of clas- 
sical contributions remains Gaussian, or perhaps Poisso- 
nian p5[ | , but the distribution for transition strengths ex- 
pected from random matrix theory is exponential (for the 
unitary ensemble) or Porter-Thomas (for the orthogonal 
ensemble). Further calculations indeed show that while 
the first and second moments agree, the higher moments 
and the full distributions differ |^ . 

Among the consequences that seem worthwhile to pur- 
sue are the dependence on the initial state and the possi- 
bility to highlight the non-Lorenzian part, the modifica- 
tion to allow for non-exponential classical escape [ p^j25| 
and the singular contributions from periodic orbits near 
bifurcations |2g|l. 
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FIG. 1. Cross section correlation function for different 
widths of the observable, increasing from M = 1 (top curve) 
to 3, 10 and 50 (bottom curve). The width of the resonances 
was uniform throughout the quantum spectrum: F = A (see 
Eq. (^s])). The continuous curves are fits to the functional 
form of Eq. (^^ with f and the coefficient of the Lorentzian 
as fitting parameters with the condition C(0) = 1. 



FIG. 2. Ratio between variance and average squared for 
observables of different widths. The different symbols corre- 
spond to different sizes of the matrix and thus to different 
values of h. 
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